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Effects of allelic variations in starch synthesis-related genes on grain quality traits 
of Korean nonglutinous rice varieties under different temperature conditions 
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National Institute of Crop Science, Rural Development Administration, Suwon 441-857, Republic of Korea 
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Influences of allelic variations in starch synthesis-related genes (SSRGs) on rice grain quality were examined. 
A total of 187 nonglutinous Korean rice varieties, consisting of 170 Japonica and 17 Tongil-type varieties, 
were grown in the field and in two greenhouse conditions. The percentages of head rice and chalky grains, 
amylose content, alkali digestion value, and rapid visco-analysis characteristics were evaluated in the three 
different environments. Among the 10 previously reported SSRG markers used in this study, seven were poly- 
morphic, and four of those showed subspecies-specific allele distributions. Six out of the seven polymorphic 
SSRG markers were significantly associated with at least one grain quality trait (R 2 > 0.1) across the three dif- 
ferent environments. However, the association level and significance were markedly lower when the analysis 
was repeated using only the 170 Japonica varieties. Similarly, the significant associations between SSRG al- 
lelic variations and changes in grain quality traits under increased temperature were largely attributable to the 
biased allele frequency between the two subpopulations. Our results suggest that within Korean Japonica vari- 
eties, these 10 major SSRG loci have been highly fixed during breeding history and variations in grain quality 
traits might be influenced by other genetic factors. 
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Introduction 

The increase in temperature as a result of global warming is 
threatening stable crop production in many countries. The 
grain yield and quality are degenerating in rice, a staple for 
half of the world's population. The global average tempera- 
ture rose by 0.85°C during 1880-2012 and an additional in- 
crease of 1.0°C-3.7°C is expected for 2081-2100 relative to 
1986-2005 (IPCC 2013). At the flowering stage, excessive 
heat causes spikelet sterility and yield reduction in rice. 
When panicles are exposed to temperature higher than 35°C 
for more than 1 hour, spikelet fertility is significantly re- 
duced in most rice varieties because of deterioration in pol- 
len viability and anther dehiscence (Jagadish et al. 2007, 
2010, Prasad et al. 2006, Rang et al. 2011, Ye et al. 2012). 
During the grain-filling stage, temperature higher than the 
optimum causes lower grain weights and decreases the per- 
centage of ripened grains because of accelerated respiration 
and delayed dry matter production and endosperm cell en- 
largement (Mohammed and Tarpley 2009, Morita et al. 
2005, She et al. 2010, Zakaria et al. 2002). Furthermore, the 
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appearance and palatability of rice are diminished under 
high temperature because of loose and irregular starch gran- 
ule formation (Chun et al. 2009, Ishimaru et al. 2009, 
Tanaka et al. 2009). Starch biosynthesis is hindered under 
high temperature, as the expressions of starch-synthesizing 
genes are suppressed, while those of starch-degrading genes 
are enhanced (Yamakawa et al. 2007). Therefore, maintain- 
ing stable yield and quality under high temperature is one 
of the most important goals in rice breeding to cope with 
climate change. 

Starch is the main component of rice, constituting 90% 
of the milled grain. Thus, the physicochemical properties of 
starch, such as amylose content (AC), gel consistency, ge- 
latinization temperature (GT), and rapid visco-analysis 
(RVA) characteristics, greatly affect the appearance and the 
processing, cooking, and eating qualities of rice grains (Bao 
et al. 2006c, Cruz et al. 1989, Fitzgerald et al. 2009, Tian 
et al. 2009). Starch biosynthesis is mainly catalyzed by the 
enzymes ADP-glucose pyrophosphorylase (AGPase), 
granule-bound starch synthase (GBSS), starch synthase (SS), 
starch branching enzyme (SBE), and starch debranching en- 
zyme (SDBE) (Jeon et al. 2010, Tian et al. 2009). Naturally, 
allelic variations in starch synthesis-related genes (SSRGs) 
affect phenotypic variations in grain quality traits. For ex- 
ample, the G/T single nucleotide polymorphism (SNP) in 
the 5' leader intron of GBSSI, which influences the splicing 
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efficiency of the gene, explains AC variations in rice along 
with the (CT) n microsatellite in the first exon of the same 
gene (Ayres et al. 1997, Bao et al. 2006a, Bligh et al. 1995, 
Dobo et al. 2010, Larkin and Park 2003). Also, the TT/GC 
SNPs in SSIIa that result in a leucine-to-phenylalanine 
substitution affect the GT and alkali digestion value (ADV) 
of rice (Bao et al. 2006b, Nakamura et al. 2005, Umemoto 
et al. 2004, Umemoto and Aoki 2005). Allelic effects of 
SBE and SDBE on grain quality traits have been document- 
ed as well (Bao et al. 2006a, Han et al. 2004, Larkin et al. 
2003). 

Recently, advanced sequencing technology has allowed 
large-scale identification of sequence polymorphisms in 
SSRGs, and their associations with the variations in grain 
quality traits have been widely investigated. By sequencing 
the entire loci of 18 SSRGs from 16 core entries of diverse 
rice germplasms and analyzing the allelic effects of SSRGs 
on starch physicochemical properties, Tian et al. (2009) iden- 
tified the roles of different SSRGs affecting rice cooking 
and eating qualities. Lee et al. (2009) sequenced complete 
or partial genie regions of eight SSRGs from 27 rice varieties 
with diverse genetic backgrounds and developed 25 cleaved 
amplified polymorphic sequence (CAPs) or derived CAPs 
markers, among which SNPs in GBSSI and AGPase-S were 
strongly associated with AC. Yan et al. (201 1) developed 43 
sequence tagged site (STS) and CAPs markers for geno- 
typing glutinous rice germplasms by sequencing 17 SSRGs 
from 13 representative varieties and studied the roles of 
SSRGs in controlling the RVA characteristics of 118 gluti- 
nous rice varieties. Through massive parallel sequencing, 
Kharabian-Masouleh et al. (2011) found 501 SNPs and 113 
insertion/deletion (InDels) markers in 17 SSRGs from 233 
Australian rice breeding lines. Among those, 66 functional 
SNPs were chosen for association analysis; their results re- 
affirmed that the GBSSI G/T SNP and SSIIa TT/GC SNPs 
are the major SSRG polymorphisms explaining variations 
in AC and GT, respectively (Kharabian-Masouleh et al. 
2012). 

While extensive studies have been conducted on the ef- 
fects of SSRG allelic variations on grain quality traits, few 
have explored the relationships between SSRG allelic varia- 
tions and the changes in grain quality traits under increased 
temperature. One of the few examples is the study of Inukai 
and Hirayama (2010), which showed that the reduction in 
AC under high temperature was greater in the genotypes 
with the Wx b allele of GBSSI than in those with the Wx* al- 
lele. Hence, more research is required to identify SSRG al- 
leles that are useful for maintaining stable grain quality un- 
der high temperature. In this study, we investigated the 
effects of allelic variations in 10 previously reported SSRG 
markers on the phenotypic variations in grain quality traits 
of 187 nonglutinous Korean rice varieties under ambient 
(field) and two greenhouse conditions. Further analyses 
were carried out to identify relationships between SSRG al- 
lelic variations and the changes in grain quality traits under 
different temperatures. 



Materials and Methods 

Plant materials and growth conditions 

A total of 187 nonglutinous Korean rice varieties consist- 
ing of 170 Japonica and 17 Tongil-type (derived from an 
Indica-Japonica cross) were used in this study. These varie- 
ties have been developed since 1979, when 'Taebaeg' was 
released, while most of the other varieties have been re- 
leased to farmers since the late 1980's. The seeds were sown 
in a seedling nursery at the Rice Experimental Station of the 
Department of Rice and Winter Cereal Crop (35.9343°N, 
126.9259°E), National Institute of Crop Science, Iksan, 
South Korea, on May 4, 2010. On June 7, 20 plants of each 
variety were transplanted in a row in the field and two 
greenhouses with a planting density of 30 cm x 15 cm, one 
plant per hill. The greenhouses were covered with the 
'F-clean® Clear film' (60 urn thickness with 94% light 
transmission; AGC GREEN-TECH Co., LTD., Tokyo, 
Japan) to mitigate sun light shadiness by covering film. The 
greenhouses were also equipped with temperature sensor- 
controlled ventilation fans at the top and manual side win- 
dows to avoid over-heating. To create different temperature 
conditions in the two greenhouses, the ventilation fans were 
set to operate when the temperature reached 32°C in green- 
house I and 34°C in greenhouse II. Side windows were 
opened to two-thirds of the full height in greenhouse I and 
half of the full height in greenhouse II. Temperatures in the 
field and the two greenhouses were recorded every 30 min- 
utes from transplanting to harvest with SK-L200T II data 
loggers (SATO, Tokyo, Japan). N, P 2 0 5 , and K 2 0 fertilizers 
were applied at 90 kg/ha, 45 kg/ha, and 57 kg/ha, respec- 
tively. The fertilization split ratio was as follows: 50% as 
basal fertilizer, 20% at the tillering stage, and 30% at the 
panicle initiation stage for N, 100% as basal fertilizer for 
P2O5, and 70% as basal fertilizer and 30% at the panicle ini- 
tiation stage for K 2 0. 

Molecular marker analysis 

Genomic DNA of the 187 rice varieties was extracted 
from the fresh young leaves of 2-week-old seedlings using 
a BioSprint 96 DNA Plant Kit (Qiagen GmbH, Hilden, 
Germany) following the manufacturer's instructions. PCR 
was carried out using AccuPower PCR Premix (Bioneer, 
Daejeon, Korea) with 10 pmol forward and reverse primers 
and 50-100 ng template DNA in a 20-fj.l total reaction vol- 
ume. The primer sequences of 10 previously reported SSRG 
markers are listed in Table 1. In a MyGenie 96 thermal 
block (Bioneer), amplifications were carried out with the 
following conditions: initial denaturation at 95°C for 3 min 
followed by 35 cycles of 45 s at 94°C, 45 s at 55-60°C, and 
1 min at 72°C, and a final extension at 72°C for 5 min. 
Banding patterns were visualized by electrophoretic sepa- 
ration of the PCR products on 1.0% agarose (AGPase-L 
InDel, SSIIa TT/GC SNPs, and SBEI InDel) or 8.0% 
polyacrylamide gels (GBSSI InDel and SDBE InDel). For 
other SSRG markers, PCR product was digested with an 
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Table 1. Starch synthesis-related gene (SSRG) markers used for genotyping 187 nonglutinous Korean rice varieties 



Gene 


Chromo- 
some 


Marker 


Sequence (5'— >3') c 


RE' 


PCR product size 
(bp) 


Reference 


AGPase-L a 


3 


InDel (TE 6 ) 


F: GCTCCATTGAGTTGAATTGGC 
R: CCGCATGGACTACATGGAC 




In: 1104 
Del: 704 


Mo era/. 2012b 


ALrrClSe-o 


o 
o 


CMP f/~7T\ 


P- ATP.TA A ArTr,TT4r,4 ATrHA A A 
r . Al VJ 1AAAL lul 1AUAA1 tuAAuA 

R: CAAGATCTATGTGCTGACACAAT 


HCOK V 


Vs. jUj 

T: 281+24 


.Lee €t al. zuuy 






CMP ( a rr\ 


p- xr.Tr.Tr.TTTr.Tr.r.TA r,r.r.Tr. a a r,r 

r. lUlUIUl 1 1U1 VJvJ 1 / W.I V.I VI 1 V.I ,' V t \ v IV 

R: CGATGCAAATTGGATTTGGATAG 


/ilu 1 


T: 315 + 23 


Lcc cl til, Z.\)\jy 


GBSSI(Wx) 


6 


SNP (T/G) 


F: CTTTGTCTATCTCAAGACAC 
R: TTTCCAGCCCAACACCTTAC 


AccI 


T: 253 

G: 127+126 


Ayres era/. 1997, 
Bao et al. 2006a 






InDel 


F: TGCAGAGATCTTCCACAGCA 
R: GCTGGTCGTCACGCTGAG 




In: 196 
Del: 173 


Wancnana et al. 2003 






CMP (T/C\ 


F- AGA AGGGGTGAGGrTTTGA ACCCGG 

r . \V II II ) I IV IV 1 1 V ! I V IV IV 111 1 1 1 . II V V VIII 

R: CTGGATGAGTCCACAGGGCTCGAAG 


Nhe I 


T: 295 

C: 244 + 51 


LCC cl HI. £.\j\jy 


SSlla 


g 




F1- rTGGATCACTTCA AGCTGTAPGAC 

J. 1 . V. 1 VJ VJ / \ 1 V I \ V 1 1 V / \ / W.I V 1 V.l 1 / \ v.. V.i t \ V 

F2: CAAGGAGAGCTGGAGGGGGC 
Rl: GCCGGCCGTGCAGATCTTAAC 
R2: ACATGCCGCGCACCTGGAAA 




TT- 1 336 + 831 
GC: 1336 + 542 


Rao pt nl ?00fih 


SBEI 


6 


InDel (TE 6 ) 


F: GAGTTGAGTTGCGTCAGATC 
R: AATGAGGTTGCTTGCTGCTG 




In: 882 
Del: 547 


Bao et al. 2006a, 
Han et al. 2004 


SBEIII 


2 


SNP (C/G) 


F: GTCTTGGACTCAGATGCTGGACTC 
R: ATGTATAACTGGCAGTTCGAACGG 


Spel 


C: 510 

G: 295 + 215 


Bao et al. 2006a, 
Han et al. 2004 


SDBE 


4 


InDel 


F: TGGCCTATCGATTGATGAGA 

R: AAAATTTTGCTGGTAATAGAAAAA 




In: 299 
Del: 280 


Larkin et al. 2003 



" L and S, large and small subunits, respectively. 

b Insertion and deletion polymorphism of a transposable element. 

c F and R, forward and reverse primers, respectively. 

d Restriction enzymes used for PCR product digestion. 

" TT/GC SNPs are two contiguous SNPs in SSIIa (Bao et al. 2006b). 



appropriate restriction enzyme before the electrophoresis on 
1.0% agarose gel (SBEIII C/G SNP) or 8.0% polyacryla- 
mide gel (AGPase-S C/T and A/T SNPs, and GBSSITIG and 
T/C SNPs). Detailed PCR and electrophoresis procedures 
for each marker are described in the references provided in 
Table 1. Based on the detected banding patterns of SSRG 
markers, allele types were assigned to each rice variety. 

Grain quality evaluation 

The 187 varieties were divided into six maturity groups 
according to their heading dates, and were harvested 112- 
129 days after transplanting, with 3-4 day intervals between 
adjacent maturity groups (Table 2). For each variety, 10 
plants with relatively uniform growth were harvested in the 
field and from the two greenhouses. The panicles were 
hand-threshed, dehulled with SY88-TH (Ssangyong, 
Incheon, Korea), and milled with VP-32T (Yamamoto, 
Yamagata, Japan) to investigate head rice percentage 
(HRP), chalky grains percentage (CGP), AC, ADV, and 
RVA characteristics. HRP and CGP were examined by scan- 
ning milled rice grain with an RN-300 Grain Inspector 
(Kett, Tokyo, Japan). AC was estimated indirectly by scan- 
ning the milled rice samples with near infrared spectroscopy 
using an Infratec 1241 Grain Analyzer (Foss, Hillerod, 
Denmark). ADV was measured visually by soaking the 
milled rice grains in 1.4% KOH solution for 24 h at 30°C as 
described by Little et al. (1958). Rice flour (3.0 g) stirred in 



Table 2. Average daily temperature during grain filling (from heading 
to harvest) in the field and two greenhouses for each maturity group 

Average daily temperature during grain-filling 



. . Field Greenhouse I Greenhouse II 

group varieties 

Min Mean Max Min Mean Max Min Mean Max 



I 


36 


21.5 


25.7 


33.4 


22.3 


26.1 


32.5 


22.8 


26.5 


33.0 


11 


27 


20.6 


24.9 


32.7 


21.5 


25.3 


31.8 


22.1 


25.9 


32.3 


III 


26 


19.5 


23.9 


31.6 


20.5 


24.4 


30.8 


21.1 


25.0 


31.3 


IV 


25 


18.8 


23.4 


31.2 


19.8 


23.9 


30.5 


20.4 


24.5 


30.8 


V 


46 


18.1 


22.8 


30.8 


19.1 


23.3 


30.2 


19.8 


23.9 


30.6 


VI 


27 


17.5 


22.3 


30.4 


18.7 


22.9 


29.9 


19.3 


23.5 


30.2 



" Plants were harvested 112-129 days after transplanting according to 
their maturity, with 3- to 4-day intervals between adjacent groups. 

25 ml water was used to measure RVA characteristics in- 
cluding pasting temperature (PT), peak viscosity (PKV), hot 
paste viscosity (HPV), cool paste viscosity (CPV), setback 
viscosity (SBV), and breakdown viscosity (BDV) with a 
Rapid Visco Analyzer (Newport Scientific Pty, Warriewood, 
NSW, Australia). 

Statistical analyses 

Descriptive statistics were calculated for grain quality 
traits under each temperature condition. The changes in 
grain quality traits under increased temperature were calcu- 
lated by subtracting the values obtained in the field from 
those obtained in greenhouse II, which had the highest 
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temperature. The percentage of phenotypic variation ex- 
plained (R 2 ) was obtained to estimate the relative contribu- 
tion of SSRG loci to grain quality traits by single-locus 
ANOVA. For F-tests, a locus with a P- value less than 0.05 
was considered to be significant. Then the least significant 
difference (LSD) test was conducted to compare the means 
of the grain quality traits according to the different allele 
types for each marker. Calculations and analyses were per- 
formed using PROC MEANS and PROC GLM of the SAS 
statistical software package (version 9.1, SAS Institute, NC, 
USA). Nei's (1973) gene diversity value for each SSRG 
marker was calculated using PopGen version 1.31 (http:// 
www.ualberta.ca/~fyeh/). To identify pair-wise genetic rela- 
tionships between rice varieties according to SSRG alleles, 
a genetic similarity matrix was calculated using simple 
matching coefficients. The similarity matrix was used to 
construct a dendrogram by the unweighted pair-group meth- 
od with the arithmetic mean (UPGMA) algorithm and the 
sequential, hierarchical, and nested clustering (SHAN) rou- 
tine in the NTSYS program (version 2.0; Rohlf 1998). 

Results 

Changes in grain quality traits under increased temperature 
The daily mean temperature from transplanting to har- 
vest in the field and two greenhouses is shown in Table 2. 
Compared to the field, the average daily mean temperature 
during grain filling was 0.4-0.6°C higher in greenhouse I 
and 0.8-1.2°C higher in greenhouse II. Also, the average 
daily mean temperature during grain filling was higher for 
early maturity groups than for late maturity groups. For ex- 
ample, the average daily mean temperature during grain 



filling in the field was 25.7°C for Group I and 22.3°C for 
Group VI. Similarly, the average daily mean temperature 
during grain filling in greenhouse II was 26.5°C for Group I 
and23.5°Cfor Group VI. 

The changes in the average values of 10 grain quality 
traits under increasing temperature in 187 nonglutinous 
Korean rice varieties are shown in Fig. 1 (see Supplemen- 
tal Table 1 for descriptive statistics). The average HRP sig- 
nificantly decreased under high temperature; it was 79.4% 
in the field, 60.2% in greenhouse I, and 49.1% in green- 
house II (Fig. 1A). There were significant increases under 
high temperature for the average CGP, which was 13.4%, 
30.7% and 39.4% in the field, greenhouse I, and greenhouse 
II, respectively (Fig. IB). Average AC showed a significant 
decrease with increasing temperature, decreasing from 
18.2% in the field to 15.4% and 14.4% in greenhouse I and 
II, respectively (Fig. 1C). There was no significant differ- 
ence in average ADV between the rice grown in the two 
greenhouses, but that of rice grown in the field was higher 
than those of rice grown in the greenhouses (Fig. ID). 
Among the six RVA characteristics, only average BDV 
showed significant changes under increasing temperature 
(Fig. 1E-1J). The average BDV was 90.8 in the field and 
significantly higher in greenhouse I (96.1) and greenhouse 
II (102.0) (Fig. 1J). The PT was significantly higher in the 
greenhouses than in the field, but there was no significant 
difference between the two greenhouses (Fig. IE). Average 
HPV, CPV, and SBV were higher in the field than in the 
greenhouses (Fig. 1G-1I). Average PKV did not change 
significantly under high temperature (Fig. IF). 

Temperature during grain filling was higher in early- 
maturing group compared to late-maturing group (Table 2). 



A 

100 

so 

60 
40 
20 
0 



Head rice (%) 



B Chalky grains (%) 

75 



C Amylosc content (%) 

25 



D Alkali digestion value 
7.5 



E Pasting temperature (°C) 
100 



F GI on 




F Peak viscosity (RVU) 

300 



240 
1X0 
120 
60 
0 



rih 



GI 



(ill 



G Hot paste viscosity (RVU) 
200 

160 

120 

80 

40 
0 



F GI Gil 



H Cool paste viscosity (RVU) I Setback viscosity (RVU) 



J Breakdown viscosity (RVU) 




Fig. 1. Changes in the grain quality traits of 187 nonglutinous Korean rice varieties under increased temperature. (A) head rice percentage 
(HRP); (B) chalky grains percentage (CGP); (C) amylose content (AC); (D) alkali digestion value (ADV); (E) pasting temperature (PT); (F) peak 
viscosity (PKV); (G) hot paste viscosity (HPV); (H) cool paste viscosity (CPV); (I) setback viscosity (SBV); (J) breakdown viscosity (BDV). 
Traits showing significant (P < 0.01) changes under increased temperature both between F&G1 and GI&GII are indicated with gray bars, while 
the others are displayed with white bars. Error bars show standard deviation. Different letters above error bars indicate significant difference at 
P < 0.01. F, GI, and Gil represent the field, greenhouse I, and greenhouse II, respectively. 
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Table 3. Effect (R 2 ) of maturity groups on the grain quality traits of 187 nonglutinous Korean rice varieties under different temperature condi- 



tions" 


11 all 




All (n = 


187)'' 






Japonica (n - 


170) c 




F 


GI 


Gil 


A 


F 


Gl 


Gil 


A 


HDD 

H K 1 


U.t /g' 1 "^ 


n ^oa* ** 


n *** 


c\ a <n*** 


U.zio 




U.Jol 


f\ A Q< * ** 


CGP 


0 142*** 


0.386*** 


0 491 *** 


0.413*** 


0.170*** 


0.427*** 


0.575*** 


0.464*** 


AC 


0.440*** 


0 474*** 


0.515*** 


0.319*** 


0.436*** 


0.478*** 


0.517*** 


0.329*** 


ADV 


Q 211*** 


0.182*** 


0.107*** 


0.118*** 


0.187*** 


0.105** 


0.067* 


0.140*** 


PT 


NS 


NS 


0.078* 


NS 


NS 


NS 


0.087* 


NS 


PKV 


0.210*** 


0.094** 


NS 


0.204*** 


0.318*** 


0.151*** 


NS 


0.206*** 


HPV 


0.354*** 


0.230*** 


NS 


0.306*** 


0.349*** 


0.239*** 


NS 


0.297*** 


CPV 


0 241*** 


0.148*** 


0.095** 


0 291*** 


0.252*** 


0.148*** 


0.092** 


0.289*** 


SBV 


NS 


0.083** 


NS 


NS 


NS 


0 129*** 


0.081* 


NS 


BDV 


0.059* 


0.067* 


NS 


NS 


0.069* 


0.076* 


0.074* 


NS 



" R 2 was calculated from ANOVA results in which the maturity group (Group I-Group VI; see Table 2) was used as an independent variable. R 2 

values are significant at *P < 0.05, **P < 0.01, and ***P < 0.001; NS no significance. 
b AC, amylose content; ADV, alkali digestion value; BDV, breakdown viscosity; CGP, chalky grains percentage; CPV, cool paste viscosity; HPV, 

hot paste viscosity; HRP, head rice percentage; PKV, peak viscosity; PT, pasting temperature; SBV, setback viscosity. 
c F, Gl, and Gil represent the field, greenhouse I, and greenhouse II, respectively. A changes in grain quality traits under increased temperature 

(calculated by subtracting the field values from their corresponding values in greenhouse 11). 



Therefore, we conducted analysis of variance (ANOVA) 
with the maturity group (I- VI) as an independent variable to 
estimate the effect of maturity on the variations in grain 
quality traits under different temperature conditions (Table 
3). There were significant associations between the maturity 
group and the variations in HRP, CGP, AC, ADV, and CPV, 
in the field, greenhouse I, and greenhouse II. The maturity 
group explained >49% of the variations in HRP, CGP, and 
AC in greenhouse II. In greenhouse II, there were notable 
differences in the three grain quality traits among maturity 
groups I-IV, V, and VI (Supplemental Table 2). The average 
HRP in greenhouse II was 37.6%-43.0% in group I-IV, 
59.2% in group V, and 71.8% in group VI, while the average 
CGP in greenhouse II was 42.8%-51.0% in group I-IV, 
30.2% in group V, and 17.2% in group VI. Similarly, the 
average AC in greenhouse II was 12.5%— 13.5% in group 
I-IV, 15.9% in group V, and 17.6% in group VI. The maturi- 
ty group also explained >31% variations in the changes in 
HRP, CGP, and AC under increased temperature {A in Table 
3). The effect of the maturity group on grain quality traits 
was also detected when the ANOVA was conducted with 
170 Japonica varieties, after excluding 17 Tongil-type 
varieties (Table 3). From this result, we could infer that the 
deterioration in grain appearance and the decrease in AC 
under high temperature were more severe in the varieties 
with early heading dates because they were exposed to 
higher temperature during the grain filling period, compared 
with those with late heading dates (Table 2). 

Allele distribution of 10 SSRG markers in 187 nonglutinous 
Korean rice varieties 

Although the maturity group (based on heading date) ex- 
plained large portion of the variations in grain quality traits 
under different temperature conditions (Table 3), there were 
also high level of variations in grain quality traits among 
different varieties within the same maturity group (Supple- 



mental Table 2). To examine the effect of SSRG alleles on 
the variations in grain quality traits under different tempera- 
ture conditions, we genotyped 187 nonglutinous Korean 
rice varieties with 10 previously reported SSRG markers 
(Table 1). 

Among the 10 SSRG markers, seven were polymorphic 
(Table 4). The SSIIa TT/GC SNPs showed the highest gene 
diversity (0.499) among the 10 SSRG markers. Although all 
Tongil-type varieties contained the TT allele, the Japonica 
varieties showed even allele distribution with 80 varieties 
carrying the TT allele and 88 varieties carrying the GC al- 
lele. AGPase-S C/T and A/T SNPs also had relatively high 
gene diversities (0.402 for C/T SNP and 0.415 for A/T SNP) 
with relatively even allele distribution among both Japonica 
and Tongil-type varieties. It should be noticed that two addi- 
tional allele types, unknown and absent, were identified by 
the DNA marker A/T SNP in AGPase-S. Lee et al. (2009) 
reported that, among 55 rice germplasm with different 
ecotypes {Japonica, Japonica/Indica hybrid, and Indica), 
the A allele produces a single band (338 bp) while the T al- 
lele cleaves into two fragments (3 15 bp and 23 bp) by appli- 
cation of an endonuclease, Alu I. However, when we geno- 
typed the Korean Japonica rice varieties, 25 varieties did 
not generate PCR products. We also observed bands with 
the size obviously smaller than 300 bp, and designated them 
as unknown (U) allele which is clearly distinguished from 
the T allele. The U allele was actually derived from a 50 bp 
deletion within the amplicon (data not shown). The results 
indicate there is sequence variations within the AGPase-S 
locus, which was not identified in the 55 rice germplasm 
evaluated by Lee et al. (2009). 

For the AGPase-L InDel, 169 out of 170 Japonica varie- 
ties had the insertion allele, while all of the 17 Tongil-type 
varieties had the deletion allele, indicating subspecies- 
specific allele distribution. This type of biased allele fre- 
quency in the two subpopulations was also observed for the 
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Table 4. Allele distribution and gene diversity of 10 SSRG markers among 187 nonglutinous Korean rice varieties according to the two subpop- 



ulations — Japonica and Tongil 



Gene 


Marker 






No. of varieties 




Nei's 


gene diversity" 




Allele 


All 


Japonica 


Tongil 


All 


Japonica 


Tongil 


AGPase-L 


InDel 


In 
Del 


169 
18 


169 
1 


0 

17 


0.174 


0.012 


0 


AGPase-S 


SNP 


C 
T 


135 
52 


127 
43 


8 
9 


0.402 


0.378 


0.498 




SNP 


A 
T 
U 6 

absent 


119 
10 

33 
25 


109 
3 
33 
25 


10 

7 
0 
0 


0.415 


0.383 


0.484 


GBSSI 


SNP 


T 
G 


187 
0 


170 
0 


17 
0 


0 


0 


0 




InDel 


In 
Del 


0 
187 


0 
170 


0 
17 


0 


0 


0 




SNP 


T 
C 


187 
0 


170 
0 


17 
0 


0 


0 


0 


SSIIa 


SNP 


TT 
GC 
absent 


97 
88 
2 


80 
88 
2 


17 
0 
0 


0.499 


0.499 


0 


SBEI 


InDel 


In 

Del 

absent 


165 
20 
2 


165 
3 
2 


0 
17 
0 


0.193 


0.035 


0 


SBEII1 


SNP 


C 
G 


18 
169 


2 
168 


16 

1 


0.174 


0.023 


0.111 


SDBE 


InDel 


In 

Del 


13 
174 


0 
170 


13 
4 


0.129 


0 


0.360 



"Nei 1973. 

b A previously unknown allele different from A and T SNPs was identified in 33 varieties and was designated as the unknown (U) allele. 



SBEI InDel, the SBEIII C/G SNP, and the SDBE InDel, all 
of which had genetic diversity values of less than 0.2. The 
three markers in GBSSI were monomorphic; all of the 187 
varieties had the T alleles for both T/G and T/C SNPs in 
GBSSI, and the deletion allele for the GBSSI InDel marker. 

A dendrogram generated from the genotyping results re- 
vealed 20 haplotypes for the 10 SSRG markers among the 
187 nonglutinous Korean rice varieties (Fig. 2). The two 
most common haplotypes were Haplotypes A and B, con- 
taining 25% and 3 1% of the tested varieties, respectively. 
The Haplotypes C, D, H, M, N, O, and T were variety- 
specific for 'Ansan', 'Cheonghaejinmi', 'Namil', 'Manjong', 
'Jinbaek', 'Nongan', and 'Dasan2', respectively. Both Haplo- 
types A and B contained varieties from all of the maturity 
groups (I- VI in Table 2). Most of the other haplotypes did 
not show maturity group-dependent variety composition ei- 
ther, indicating that the allele distributions of most SSRG 
markers tested in this study are likely to be independent 
from heading dates. In contrast, the haplotypes of the 170 
Japonica varieties (A-O) were clearly distinguished from 
those of the 17 Tongil-type varieties (P-T) in the dendro- 
gram, reflecting the subspecies-specific allele distribution of 
many SSRG markers, i.e., AGPase-L InDel, SBEI InDel, 
SBEIII C/G SNP, and SDBE InDel (Fig. 2 and Table 4). 



Effects of SSRG allelic variations on grain quality traits 
under different temperature conditions 

Table 5 shows the significant associations between the 
allelic variations in SSRG markers and the phenotypic vari- 
ations in grain quality traits of 187 nonglutinous Korean 
rice varieties under three different temperature conditions 
(the field, greenhouse I, and greenhouse II). Among the 
seven polymorphic SSRG markers, six showed a relatively 
high level of association (R 2 > 0.1) with variations in one or 
more grain quality traits of the 187 varieties in the three dif- 
ferent environments; i.e., AGPase-L InDel (PKV, SBV, 
BDV), AGPase-S A/T SNP (SBV, BDV), SSIIa TT/GC 
SNPs (ADV), SBEI InDel (SBV, BDV), SBEIII C/G SNP 
(PKV, SBV, BDV), and SDBE InDel (PKV, SBV, BDV). 

The data were reanalyzed after excluding 1 7 Tongil-type 
varieties to rule out the effect of population structure. In this 
analysis, significant associations were not observed in most 
markers, or the R 2 value decreased to a large extent 
(Table 5). For example, although allelic variations in the 
AGPase-L InDel significantly explained some of the pheno- 
typic variations in CGP (2.1-13.8%), AC (2.3-5.7%), ADV 
(3.8-6.9%), and PKV (13.1-29.1%) of the 187 varieties in 
the three different environments, there was no significant 
relationship when the 17 Tongil-type varieties were exclud- 
ed from the analysis. Also, while more than 29% of the var- 
iations in SBV and BDV of the 187 varieties was explained 
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Coefficient 
0.6 0.7 oi 



0.456 



Haplotype 
(No.) 

A (46) 



B(58) 




F(10) 

G(3) 
H,I(4) 



K(19) 



1 -Boseok 
1 -Dunnae 
l-Geumo3 
1-Hanseol 
l-Jinbuol 
I-Joun 

1-Junghwa 1-Urani 
1 -Manlio 2-Cheonga 



1 -Munjang 
1 -Namweon 
I-Odae 
I-Odael 
I-Sangju 



2-Daejin 2-Naepung 3-Kwangan 5-Dongan 

2-Hwadong 2-Obong 3-Migwang 5-Geumnam 

2-Hwangkeumbora 3-Anseong 3-Seoan 5-Namgang 

2-Jinmi 3-Cheongnam 4-Cheongmyeong 5-Weonhwang 

2-Jinpum 3-Gancheok 4-DaeripI 6-Hoan 



1-Shinunbong 2-Joryeong 
2-Jungan 
2-Jimgsan 



3-Geuman 

3-Hwabong 

3-Hwaseong 



4-Hwajung 

4-Seoanl 

4-Seojin 



6-Saenuri 



I -Geumyeong 1-Pyeongweon 2-Saesangju 3-Sampyeong 
2- S anibaeg 3 - Sura 
2-Sandeuljinmi 4-Cheongan 



1-Goun 
1-Gru 
1-Iiiweol 
1-Jinbong 
1-Jinbu 
I -Joan 
1-Kcumo 
1 -Marian 
1-Manchu 



1-Samcheon 
]-Sinunbongl 
1-Sobaeg 
]-Taebong 
1 -Taeseong 
1-Unbong 
]-Undoo 

1- Unjang 

2- Handeul 



4-Nagdong 5-Jongnam 

4- Seogjeong 5-Mangeum 

5- Cheongcheongjinmi 5-Manweol 



2-Sangini 

2- Sangsan 

3- Anjung 
3-Haeddeumi 
3-Jangan 

3- Juan 
3-Manpung 



4-Chilbo 

4-Daean 

4-Daeya 

4-Donghae 

4-Hwaan 

4-Hwajin 



5-Daecheong 

5-Dami 

5-Gopum 

5-Gyehwa 

5-Haiami 

5-Hwasbinl 



5-Mibyang 
5-Nampyeong 

5- Palgong 

6- Dongjin 
6-Hojin 



4-Hwamyeong 5-Hyangnam 



2-Ansan 6-Cheonghaejinmi 



3-Kxumo2 3-Yeonghae 5-Hwashin 6-Hopyimg 6-llmi 6-Yeongnani 



1 -Jogwang 5-Juanl 6-Daesan 6-Manmi 6-Tamjin 
5-Ilpum 6-Chucheong 6-Hwanam 6-Nongho 6-Yeongbojinmi 



2-Manna 4-Nogyang 5-Seomyeong 



2-Namil 5-Hopum 5-Junam 6-Dacheong 6-Hwarang 



1- Junamjosaeng 5-Boramcban 5-Donghaejinmi 5-Hwanggeurnnodeul 5-Seogan 

2- Unkwang 5-Borami 5-Dongjin2 5-Onnuri 5-Sindongjin 

3- Sobi 5-Clieongho 5-Hanmaeum 5-Pyimgan 



2- Pungmi 3-Haepyeong 4-Samdeog 5-Dongjin I 5-Pungmi 1 

3- Daepyeong 3-Samgwang 4-Sangok 5-Hanam 5-Saegyehwa 
3-Gangbaeg 4-Haechanmulgyeol 4-Seopyeong 5-Jinsumi 5-Yeongan 
3-Haeoreumi 4-Hwayeong 5-Deuiaechan 5-Malgemni 



j 6-Honong 6-Hwangkeumnuri 6-Hwasam 6-Sujin 
O 2-Manjong 6-Jinbaek 3-Nongan 



2-Hyangmi2 3-Hyangmil 4-Areum 4-Nampung 
2-Taebaeg 4-Samgang 6-Segyejinmi 
5-Namcheon 6-Jangseong 



2-Yongju 2-Yongmoon 4-Dasanl 5-Dasan 6-Anda 6-Hanareum 6-Keunseom 



Fig. 2. Dendrogram of 187 nonglutinous Korean rice varieties by UPGMA cluster analysis based on 10 SSRG loci. A total of 20 unique SSRG 
haplotype groups were identified within the 187 rice varieties. For each haplotype group, the number of rice varieties is shown in parentheses, 
except for haplotype groups with only one entry. Variety names are shown beside the corresponding haplotype group; Arabic numeral prefixes 
show their designated maturity groups (see Table 2). Note that the node in gray (0.456) split 187 nonglutinous Korean rice varieties into two sub- 
species — 170 Japonica varieties and 17 Tongil-type varieties, in which 15 (A-O) and 5 (P-T) SSRG haplotypes were identified, respectively. 



by the AGPase-L InDel polymorphism in the three different 
environments, less than 4% was explained after excluding 
the 17 Tongil-type varieties. Similar trends were observed 
for the AGPase-S C/T SNP (SBV, BDV), the AGPase-S All 
SNP (SBV, BDV), the SBEI InDel (ADV, PKV, SBV, BDV), 
the SBEIII C/G SNP (PKV, SBV, BDV), and the SDBE 
InDel (AC, SBV, BDV). 

The SSIIa TT/GC SNP marker was the only SSRG mark- 
er showing a relatively high level of association with certain 
grain quality traits in the field and two greenhouse condi- 
tions when the analysis included only the 170 Japonica 
varieties (Table 5). It explained 12.6%-34.2% of the ADV 
variation in 170 Japonica varieties, which was similar to 
the association analysis conducted without considering the 
population structure comprising Japonica and Tongil-type 
varieties. 

Relationship between SSRG allelic variations and changes 
in grain quality traits under increased temperature 

To evaluate the allelic effects of SSRGs on maintaining 



stable grain quality traits under high temperature, we ana- 
lyzed the relationships among SSRG allelic variations and 
the changes in grain quality traits under increased tempera- 
ture (Table 6). The changes in grain quality traits under in- 
creased temperature were calculated by subtracting the field 
values from their corresponding values in greenhouse II, 
which had the highest temperature (Table 2). 

There were significant differences in the changes in HRP 
and CGP under increased temperature (zlHRP and zlCGP, 
respectively) among the varieties with different alleles of 
six SSRGs (Table 6). zlHRP and ACG? were significantly 
lower in the varieties containing the deletion allele for 
AGPase-L InDel, the T alleles for both C/T and AT SNPs in 
AGPase-S, the deletion allele for SBEI InDel, the C allele 
for SBEIII C/G SNP, and the insertion allele for SDBE 
InDel. The AGPase-S C/T SNP was the only locus showing 
significant relationship with the change in AC under in- 
creased temperature (zlAC); A AC was significantly lower in 
the varieties carrying the T allele than in those carrying the 
C allele. The AGPase-L InDel, SBEI InDel, and SBEIII C/G 
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Table 5. Summary of association analyses (R 2 ) between SSRG marker alleles and grain quality traits among 187 nonglutinous Korean rice vari- 
eties under different temperature conditions 



Gene (Marker 3 ) 


Trait 4 




All (n = 187)' 






Japonica (11 — 170) £ 




F 


(jl 


■ i i 
ull 


F 


(jl 


(_rll 


AGPase-L (InDel) 


CGP 


0.021* 


0.078*** 


0.138*** 


NS 


NS 


NS 




AC 


O.UZJ* 


0.U5 


A AO A * 

0.U24* 


NS 


NS 


NS 




ADV 


0.045** 


0.038** 


0.069*** 


NS 


NS 


NS 




PKV 


0.189*** 


0.291*** 


0.131*** 


NS 


NS 


NS 




SBV 


0.315*** 


0.313*** 


0 297*** 


0.030* 


0.038* 


0.037* 




BDV 


0.432*** 


0.418*** 


0.344*** 


0.035* 


0.035* 


0.029* 


AGPase-S (C/T) 


SBV 


0.074*** 


0.035* 


0.034* 


0.046** 


NS 


NS 




BDV 


0.081*** 


0.104*** 


0.074*** 


0.044** 


0.078*** 


0.038* 


AGPase-S (A/T) 


SBV 


0 152*** 


0 123*** 


0 119*** 


NS 


NS 


NS 




BDV 


0.178*** 


0 144*** 


0 134*** 


NS 


0.065** 


NS 


SSJla ( 1 1/CjL) 


ADV 


U.28 /*** 


A n A A * * * 


0 1 4 1 * * * 


0.277** 


<* o 342*** 






PT 


0.110*** 


0.043** 


0.050** 


0.104*** 0.054** 


0.041** 


SBEI (InDel) 


ADV 


0.078*** 


0.049** 


0.093*** 


0.037* 


NS 


0.033* 




PKV 


0.135*** 


0.188*** 


0.070*** 


NS 


0.028* 


0.032* 




SBV 


0.215*** 


0.205*** 


0.198*** 


NS 


NS 


NS 




BUV 


A 1 1 1 * ** 

U.J 1 j^^^ 




U.zJo 


NS 


NS 




SBEIII (C/G) 


PKV 


0.156*** 


0 224*** 


0.100*** 


NS 


NS 


NS 




SBV 


0.284*** 


0.275*** 


0.261*** 


0.029* 


NS 


NS 




BDV 


0.372*** 


0.352*** 


0 292*** 


NS 


NS 


NS 


SDBE (InDel) 


AC 


0.026* 


0.063*** 


0.033* 










ADV 


0.076*** 


0.110*** 


0.110*** 










PKV 


0.157*** 


0.258*** 


0 127*** 










SBV 


0 293*** 


0 214*** 


0.213*** 










BDV 


0.312*** 


0.320*** 


0.269*** 









* Number of allele types for each marker is listed in Table 4. 

b Traits listed had significant R 2 values in 187 rice varieties under the three different conditions. Traits showing relatively high level of association 
(R 2 > 0.1) with the corresponding SSRG marker in the three different conditions are highlighted in bold. AC, amylose content; ADV, alkali di- 
gestion value; BDV, breakdown viscosity; CGP, chalky grains percentage; PKV, peak viscosity; PT, pasting temperature; SBV, setback viscosi- 
ty; R 2 values are significant at *P < 0.05, **P < 0.01, and ***P < 0.001; NS no significance; - monomorphism. 

0 F, Gl, and Gil represent the field, greenhouse I, and greenhouse II, respectively. 



SNP had significant associations with the change in BDV 
under increased temperature (zIBDV); rice varieties with 
deletion, deletion, and C alleles for the three markers had 
significantly lower zfBDV than those with the other alleles. 
When the analyses were carried out only for the Japonica 
group (n = 170), however, the significant associations dis- 
appeared or were markedly decreased for the AGPase-L 
InDel, the SBEI InDel, the SBEIII C/G SNP, and the SDBE 
InDel (Tables 6, 7). Meanwhile AGPase-S C/T SNP showed 
significant association with zfHRP, zfCGP, and AAC, and 
AGPase-S A/T SNP exhibited significant association with 
zlHRP and zlCGP in the 170 Japonica varieties, which was 
similar to the association analyzed in the total of 187 varie- 
ties. However, their effects on the grain quality traits were 
minor (R 2 < 0.07, Table 7). 

Discussion 

Among the 10 grain quality traits investigated in this study, 
four showed significant changes with increasing tempera- 
ture. Grain appearance deteriorated as the average HRP de- 
creased, while the average CGP increased under higher 
temperatures. This was similar to the previous findings that 



the deterioration in grain appearance under high tempera- 
ture was mainly because of the increase in chalkiness (Mo 
et al. 2012a). It was reported that AC decreased under high 
temperature in most rice varieties, although a few varieties 
showed increase or no significant change (Cruz et al. 1989, 
Tanaka et al. 2009, Zhong et al. 2005). Similarly, the aver- 
age AC of the 187 nonglutinous Korean rice varieties de- 
creased significantly as temperature increased. It is assumed 
that suppression of starch biosynthesis genes such as GBSSI 
and SBE under high temperature might explain the decrease 
in AC (Yamakawa et al. 2007). Among the six RVA char- 
acteristics, only BDV showed significant change under 
increased temperature. Tanaka et al. (2009) reported that 
PKV and BDV increased while SBV decreased significantly 
under high temperature in three Japonica. rice varieties. 
Zhong et al. (2005) showed that CPV and SBV decreased 
significantly under high temperature in four Indica rice vari- 
eties. In this study, the average BDV of the 187 varieties in- 
creased significantly as the temperature rose while the other 
RVA characteristics and ADV did not change significantly 
in response to higher temperature. Therefore, we decided to 
trace the changes in HRP, CGP, AC, and BDV to estimate 
the levels of stable grain quality traits under increased 
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Table 6. Changes in grain quality traits in 187 nonglutinous Korean 
rice varieties under increased temperature according to allelic differ- 
ences in SSRG markers" 



Gene 


Allele 


No. of 


zlHRP 


zlCGP 


zlAC 


zlBDV 




(0/ \ 

(/o) 


(0/ \ 

(%) 


(%J 


(KV U ) 


variety 


AGPase-L 


In 


169 


-32.2a 


27.4a 


-3.9a 


12.0a 




Del 


18 


-15.9b 


9.7b 


-3.3a 


3.8b 


AGPase-S 


C 


135 


-32.9a 


28.2a 


-4.1a 


11.6a 




T 


52 


-24.7b 


19.3b 


-3.3b 


10.0a 


AGPase-S 


A 


119 


-34.3a 


30.0a 


^1.3a 


10.8a 




T 


10 


-20.7b 


12.6b 


-3.7a 


3.7a 




U 


33 


-29.7a 


24.7a 


-3.6a 


11.9a 


SSIIa 


TT 


97 


-28.9a 


24.0a 




11.2a 




GC 


88 


-32.4a 


27.6a 


-3.8a 


10.8a 


SBEI 


In 


165 


-32.1a 


27.5a 


-3.9a 


11.8a 




Del 


20 


-17.6b 


11.2b 


-3.6a 


4.6b 


SBEIII 


C 


18 


-18.4a 


7.0a 


-3.9a 


4.1a 




G 


169 


-31.9b 


27.7b 


-3.9a 


11.9b 


SDBE 


In 


13 


-12.1a 


6.6a 


-3.1a 


4.9a 




Del 


174 


-32.0b 


27.2b 


-3.9a 


11.7a 



" Changes in grain quality traits under high temperature were calculat- 
ed by subtracting the field values from their corresponding values in 
greenhouse II. Only traits showing significant changes under in- 
creased temperature are listed (see Fig. 1). Values with a different 
letter are significantly different at P < 0.05 (highlighted in bold). 
Traits showing significant allelic effect in Japonica varieties (n = 
170) are underlined. AC, amylose content; BDV, breakdown viscosi- 
ty; CGP, chalky grains percentage; HRP, head rice percentage. 

temperature in Korean nonglutinous rice varieties. 

The early maturing varieties were exposed to higher tem- 
perature during the grain-filling period compared to the 
late-maturing varieties. Therefore, a large proportion of the 
variations in grain quality traits under different temperature 
conditions was explained by the maturity group (Table 3). 
However, large variations in grain quality traits were ob- 
served even among the varieties in the same maturity group 
(Supplemental Table 2). Therefore, we investigated the ef- 
fect of previously reported SSRG loci on grain quality traits 
under different temperature conditions to identify SSRG al- 
leles favorable for maintaining grain quality under high 
temperature. 

Allelic variations in GBSSI are well-known for determin- 
ing starch physicochemical properties such as AC and RVA 
characteristics. These characteristics affect the processing, 
cooking, and eating qualities of rice (Dobo et al. 2010, 
Kharabian-Masouleh et al. 2012, Larkin et al. 2003, Larkin 
and Park 2003, Tian et al. 2009, Yamanaka et al. 2004). 
However, there was no polymorphism for the three markers 
in GBSSI (T/G SNP, T/C SNP, and InDel) among the 187 
nonglutinous Korean rice varieties used in this study. As for 
the GBSSI T/G SNP, Bao et al. (2006a) revealed that the 
Wx a (G allele) and Wx h (T allele) are evenly distributed 
among Chinese nonglutinous rice breeding lines, reflecting 
the breeders' selections for varieties with various eating and 
cooking qualities. In contrast, all of the nonglutinous Kore- 
an rice varieties used in this study had only the Wx h allele. 
Considering that Wx a and Wx h are the predominant alleles 
for Indica and Japonica subspecies, respectively (Hirano et 



Table 7. Association (R 2 ) of SSRG allelic variations with changes in 
grain quality traits of nonglutinous Korean rice varieties under in- 
creased temperature 



Marker" 


Trait b 


All 
(n= 187) 


Japonica 
(n= 170) 


/i Kjrus tr-i> muei 


1 II P P 

Zlrllvr 








zfCGP 


0.106*** 


NS 




zlBDV 


0.029* 


NS 




/(HRP 
Zlrllvr 


n n^o*** 
yj.ujy 


ft 




zlCGP 


0.061*** 


0.043** 




ZlAC 


0.031* 


0.042** 


AGPase-S A/T 


zfHRP 


0.056* 


0.063** 




zfCGP 


0.079** 


0.052* 


SBEI InDel 


zfHRP 


0.088*** 


NS 




zlCGP 


0.098*** 


NS 




zlBDV 


0.025* 


NS 


SBEIII C/G 


zlHRP 


0.069*** 


NS 




zfCGP 


0.145*** 


0.066*** 




zfBDV 


0.027* 


NS 


SDBE InDel 


zlHRP 
zlCGP 


0.111*** 
0.106*** 





" Only the markers significantly associated with the changes in grain 
quality traits under increased temperature were analyzed. The num- 
ber of allele types for each marker is listed in Table 4. 

b AC, amylose content; BDV, breakdown viscosity; CGP, chalky 
grains percentage; HRP, head rice percentage. R 2 values are signifi- 
cant at *P < 0.05, **P < 0.01, and ***P < 0.001 ; NS no significance, 
- no polymorphism. 

al. 1998, Sano et al. 1986), it was notable that all of the 17 
Tongil-type (derived from Indica-Japonica cross) varieties 
contained only the Wx h allele. This implied that the Korean 
consumers' preference for a soft and sticky texture of 
Japonica cooked rice has been reflected in the breeders' 
selections in developing Tongil-type varieties. According to 
previous reports, GBSSI expression and AC decreased 
under high temperature to a greater extent in varieties with 
the Wx h allele than in those with the Wx a allele (Hirano 
and Sano 1998, Inukai and Hirayama 2010, Larkin and 
Park 1999). Further research is underway to introduce the 
Wx a allele into elite Korean Japonica varieties by marker- 
assisted backcrosses. This will allow us to evaluate and ex- 
ploit the effect of Wx a to retain stable grain quality under 
high temperature, to cope with global warming. 

Among the seven polymorphic SSRG markers, four 
showed subspecies-specific allele distribution in the 187 
nonglutinous Korean rice varieties. The SBEI InDel had 
significant associations with ADV, PKV, SBV, and BDV in 
the field and two greenhouse conditions and SBEIII C/G 
SNP had significant associations with PKV, SBV, and BDV 
in the three different environments. This finding was similar 
to the previous reports that the SBEI InDel and SBEIII C/G 
SNP markers were significantly associated with variations 
in RVA characteristics (Bao et al. 2006a, Han et al. 2004). 
The AGPase-L InDel and SDBE InDel were also relatively 
strongly associated (R 2 > 0.1) with PKV, SBV, and BDV in 
the field and two greenhouse conditions. However, when the 
association analysis was conducted after separating Japonica 
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and Tongil-type varieties, the significance level was highly 
reduced or disappeared because of the biased allele frequen- 
cy between the two subpopulations (Table 5). 

Only three out of 10 SSRG markers showed relatively 
even allele distributions among the Japonica varieties. The 
C/T and A/T SNPs in AGPase-S were relatively evenly dis- 
tributed among both Japonica and Tongil-type varieties 
(Table 4). Lee et al. (2009) reported that the two SNPs in 
AGPase-S explained >40% of the AC variation in 55 rice 
accessions with diverse geographical origins; however, no 
significant association with AC was found in the present 
study. Instead, the AGPase-S C/T SNP was significantly as- 
sociated with SBV and BDV and the A/T SNP was signifi- 
cantly associated with SBV and BDV in the field and two 
greenhouse conditions. Nevertheless, when the analysis was 
conducted on the 170 Japonica varieties, the significant as- 
sociations were markedly decreased. The SSIIa TT/GC 
SNPs were also evenly distributed among the Japonica va- 
rieties, with 80 varieties carrying the TT allele and 88 carry- 
ing the GC allele. This was consistent with the study by Bao 
et al. (2006b), in which 334 Chinese rice varieties also 
showed even allele distributions of SSIIa TT/GC SNPs. 
Also, the strong associations between SSIIa TT/GC SNPs 
and ADV and PT reported in other studies were verified in 
this study in all three environments (Bao et al. 2006b, 
Nakamura et al. 2005, Umemoto and Aoki 2005). In par- 
ticular, we observed a significant association (R 2 > 0.126) 
between SSIIa TT/GC SNPs and the variation in ADV in the 
Japonica group, which was similar to the association ana- 
lyzed in the total of 187 varieties. In summary, out of the 
10 SSRG markers investigated in this study, only the two 
SNPs (C/T and A/T) in AGPase-S showed relatively even 
allele distributions in both Japonica and Tongil-type non- 
glutinous Korean rice varieties. The SSIIa TT/GC SNPs 
were the only marker showing significant associations 
(R 2 > 0.1) with a grain quality trait in the three different en- 
vironments for the 170 Japonica varieties. 

In the association analysis between the SSRG allelic var- 
iations and the changes in grain quality traits under in- 
creased temperatures, all polymorphic SSRG markers ex- 
cept for SSIIa TT/GC SNPs showed a significant allelic 
effect on certain grain quality traits (Table 6). However, 
when the analysis was conducted only on Japonica varie- 
ties, only the C/T and A/T SNPs in AGPase-S showed sig- 
nificant allelic effects on the variations in certain grain 
quality traits by increased temperature, which was similar to 
the association analysis conducted without excluding the 17 
Tongil-type varieties (Table 7). Since the varieties with the 
T alleles for both C/T and A/T SNPs in AGPase-S had sig- 
nificantly less decrease in HRP and less increase in CGP 
under increased temperature, these alleles could be used in a 
breeding program to develop rice varieties that retain stable 
grain appearance under high temperature. However, it is 
likely that their effect would be limited because the two 
markers explained only 4.3-6.3% of the variations in zlHRP 
and zlCGP. As for the four subspecies-specific SSRG mark- 



ers, the changes in the grain quality traits under increased 
temperature were significantly smaller in the varieties with 
7n<i/ca-specific alleles than in those with Japonica-specific 
alleles. This is reasonable considering that Indica sub- 
species originated from the tropics while the Japonica sub- 
species are more adapted to a temperate climate. Significant 
allelic effects of the four subspecies-specific SSRG markers 
largely disappeared when the analysis was conducted only 
on the 170 Japonica varieties (Table 7). Jin et al. (2010) 
mentioned that functional loci with subspecies-specific al- 
lele distribution would not be detected if population struc- 
ture is considered in the association analysis. Therefore, fur- 
ther research is needed to confirm whether the effect of the 
subspecies-specific SSRG alleles on the changes in grain 
quality traits under increased temperature is merely because 
of the population structure or not. Overall, the level of poly- 
morphism in the 10 SSRG markers was highly limited or 
subspecies-specific among the 187 nonglutinous Korean 
rice varieties. Accordingly, the polymorphisms had limited 
abilities to explain variations in grain quality traits under 
different temperature conditions. It is assumed that the 10 
SSRG loci have been highly fixed with optimum alleles 
according to market demands during the history of Korean 
rice breeding. Therefore, there may be other genetic factors 
influencing the grain quality variations among the Korean 
Japonica varieties under increasing temperature. Another 
association study is currently being carried out with genome- 
wide SNP and SSR markers. We hope that this will allow us 
to identify other genetic factors related to the variations in 
grain quality traits among Korean nonglutinous Japonica 
rice varieties under different temperature conditions. Efforts 
to improve genetic diversity in terms of grain quality under 
high temperature are also being made using germplasm 
search and mutagenesis approaches. 
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